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We theoretically investigate the transmission of electromagnetic radiation through a metal plate 
with a zero-e metamaterial slit, where the permittivity tends towards zero over a given bandwidth. 
Our analytic results demonstrate that the transmission coefficient can be substantial for a broad 
range of slit geometries, including subwavelength widths that are many wavelengths long. This novel 
resonant effect has features quite unlike the Fabry-Perot-like resonances that have been observed 
in conductors with deep channels. We further reveal that ultranarrow zero-e channels can have 
significantly greater transmission compared to slits with no wave impedance difference across them. 



With the current state of the art in nanofabrication 
technologies, and the recent observation of resonant opti- 
cal transmission through small metal holes [l[ , there has 
been a renewed interest in electromagnetic wave trans- 
mission and diffraction in metallic nanostructures. Al- 
though diffraction effects involving optical slits and grat- 
ings has a long history, scientific curiosity coupled with 
the potential for device applications has spurred consid- 
erable research activity involving the manipulation and 
confinement of electromagnetic waves in nanoscale reso- 
nant structures. With the concurrent advent of metama- 
terials, or composite structures with tailored electromag- 
netic properties, the interplay of classical transmission 
systems with these new composite media has become an 
important issue. By engineering metallic nano-apertures, 
gratings, or channels to incorporate metamaterials, it is 
anticipated that various enhanced or exotic transmission 
characteristics will ensue. 

Lately, there has been interest in zero-e metamaterials, 
which are structures that exhibit an effective permittivity 
of zero (or nearly so) in the passband, akin to the prop- 
erty of some precious metals Q near their plasma fre- 
quency. It was shown that waveguide devices containing 
zero-e inclusions could potentially be more efficient via 
the reduction of unwanted reflections if one of the phys- 
ical dimensions was made smaller Q. This can possibly 
negate the ill effects of impedance mismatch. A multi- 
layer structure having nonlinear electric and magnetic re- 
sponses, and a refractive index close to zero, was shown 
to effectively shield [3] electromagnetic fields. Experi- 
mental work has verified @ electromagnetic tunneling 
through zero-e metamaterials at microwave frequencies. 
It was shown |7| that the effective magnetic permeabil- 
ity, fj,, can also be resonantly tuned to vanish, creating 
a matched metamaterial with an effective zero index of 
refraction. These media have been linked to applications 
in miniaturized resonators, highly directive antennas 
delay lines with zero-phase difference I/O, and transform- 
ers that convert small-curvature wave fronts into output 
beams with planar-like wave fronts 0. 

A fundamental system in which to investigate zero-e 
diffraction and transmission effects is a subwavelength 
channel through a metal film. If the slit is filled with air, 
it has been established that for light polarized perpen- 



dicular to the slit in an optically thick perfect metal of 
length, I, Fabry-Perot-like waveguide modes arise when 
I is approximately a half integer number of wavelengths: 
I w nX/2. These harmonic modes, which follow from ge- 
ometrical arguments, result in transmission peaks when 
the waves coherently superimpose over the given path 
length. As the geometrical parameters vary, the resonant 
wavelength can shift [13, EJ in metals with slit perfora- 
tions, as observed with microwaves [T3|. For ultrashort 
incident pulses, the Fabry-Perot-like modes can be res- 
onantly activated (l3l | through the Fourier components 
of the wave packet inside the slit, leading in some cases 
to enhanced transmission. The resonant enhancement of 
a nanometer scale electromagnetic pulse was also shown 
to be spatially and temporally localized [3] in the near 
field. For perfect metals with holes rather than slits, 
transmission is strictly limited to incident wavelengths 
that are less than twice the diameter of the openings. 

In this Letter we reveal some unexpected and exotic 
transmission behavior of light through a subwavelength 
zero-e slit. Our theoretical framework demonstrates that 
significant transmission can occur in these structures for 
a considerable range of ultranarrow widths that are suf- 
ficiently deep. The inverse relationship between the slit 
length and width that achieves maximal transmission is 
shown to be highly nontrivial. We show that zero-e meta- 
materials with inherently large intrinsic impedances, can 
have greater transmission than matched zero index slits, 
where n is also vanishingly small. We also take advantage 
of the subwavelength geometry (that is, the slit width 
is smaller than about half the incident wavelength), to 
invoke the single mode approximation, which has been 
shown to yield valuable physical insight into resonance 
phenomena for perfect [15[ and real [16| metals. Our an- 
alytic results thus permit one to efficiently map out the 
relevant sector of parameter space deemed appropriate 
for any future experimental endeavors. 

We consider a planar perfect metallic structure that is 
translationally invariant in the x — z plane and has chan- 
nel length I normal to the plane with width p (along a;). 
The incident beam is TM polarized, so that the magnetic 
field H is directed along z. The wavevector, fe, forms an 
angle 6 with the normal to the plane. For this configu- 
ration, the z component of the magnetic field, H z , must 
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satisfy the scalar Helmholtz equation 
8 2 H Z d 2 H 



dx 2 



Q y 2 ' 6 V k H z = 0, 



(1) 



where fco = u>/c. Our focus is primarily metamate- 
rial slits very near the characteristic plasma frequency 
and, as in recent works, represent the frequency dis- 
persive electrical response by an effective Drude model, 
e = l — Wp/[w(u + ir)], where T is correlated to the mean 
free path in the filling material. Other than when dis- 
cussing absorption loss, we generally take T equal to zero 
in order to isolate the electromagnetic effects inherent to 
vanishing constitutive relations. 

We solve ([T|) in a given region, and then match each 
solution at the corresponding interfaces using the appro- 
priate boundary conditions. To construct the solution in 
the "continuum" regions of free space surrounding the 
metal sheet, we Fourier transform Eq. (fTJ) along x and 
use separation of variables. This renders an expression 
written in terms of plane wave expansions, 
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where kn x = ko sin 9, ko y = ko cos 9, and (3 = \/k 2 — a 2 . 
Here the subscripts 1 and 3 denote the entrance and exit 
regions respectively. The unknown coefficients R(a) and 
T(a) are determined below. The general eigenmode ex- 
pansion for the field within the slit region is [Io| , 



H z , 2 = ™s[Mz - P/2)] [a m e lTmV + b m e~ lT ^\ , 



m=Q 



(4) 



where k m = mir/p, and r m = -^efik 2 — k^. The summa- 
tion is over a single index m reflecting the interconnec- 
tion among the wavenumbers consistent with the wave 
equation. The electric field is obtained directly via ap- 
plication of Maxwell's equation: E = i/(efco)V x H. It 
is generally valid to retain only the lowest order mode in 
the expansion (j4|) for the ultranarrow channels considered 
in this paper. Moreover, as e goes to zero, the summa- 
tion results in higher order evanescent fields that have a 
negligible contribution to the time-averaged energy flux. 
Thus, for the electric field we have, 
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where r = r , rj — y/fi/e is the intrinsic impedance, 
and the function £(x) is unity in the slit and vanishes 
elsewhere. 

Matching the electric field at the entrance and exit 
interfaces, yields the R and T coefficients: 
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FIG. 1: Real and imaginary components of Io as a function 
of the dimensionless width q = vrp/X. The approximate solu- 
tions and exact integral [Eq. |(HJ] have good agreement over 
the relevant range of widths. As q vanishes, the imaginary 
component has a slow divergence [see Eq. |(9j]. 



where we define $(x) = \J 2 / it sm.{x / 2) / x . A convenient 
technique [l5| that utilizes the boundary conditions to 
determine the remaining unknown a and b coefficients 
is the use of the following relationships at the appropri- 
ate openings: {H Zt i(x,0),((x)) = (H Zj2 (x,Q),((x)), and 
(H Zl2 (x,-l),((x)) = (H Zj3 (x,-l),C(x)}. After some te- 
dious algebra, we have for b, 
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where w = koprj, £ = rl, and a = be 21 ^ (wIq — 1) / (wIq+1) ■ 
The complex valued integral, Iq is a function solely of q: 
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where q = knp/2 is the dimensionless measure character- 
izing the slit width. The series expansion for the imagi- 
nary component of Iq has the following leading terms: 
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where 7 is Euler's constant. Similarly for the real part 
of I , we have (to fourth order), Re [Jo] ~ \ — j^q 2 + 
g 4 . It will be seen shortly that Iq is a very important 
quantity in the determination of the transmission. The 
accuracy of these expansions is demonstrated graphically 
in Fig. [1] where the truncated expansions are plotted 
alongside the numerically integrated Jo [Eq. ([8])]. Clearly, 
the approximations for In are satisfactory for the small q 
of interest here, deviating only slightly for the larger q, 
as would be expected for a power series centered about 
the origin. 

Next, we insert the calculated a and b coefficients into 
Eq. (O, to determine the electric field at the bottom of 
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the slit, 
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which is valid for arbitrary /j, and e. The magnetic field 
is easily obtained via the relationship, = 2qIoE x ot . 

The transmission coefficient, r, is defined as the time- 
averaged Poynting vector at the bottom of the slit, 
(Sy) y =-i, integrated over the exit opening, and divided 
by the Poynting vector of the incident plane wave: 
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which is dimcnsionless after incorporating a normaliza- 
tion factor, fco- The angular dependence is encapsulated 
entirely by the expression, T{6) = sine 2 (q sin 0) / cos 0, 
which depends weakly on q for subwavelength slits, and 
is close to unity for source waves near normal inci- 
dence. For small q, we can further write (for sin£ ^ 
0): t w 8(/j,/e)q 2 csc 2 (koly/ejl), otherwise if we consider 
channels with some integer multiple of the Fabry-Perot 
length, I = \/{2yfeJL), Eq. (fTTj) approximately reduces 
to, t w i j Jq [ 2 , which clearly only depends on the ratio 
of the slit width to the wavelength [see Eq. J§J|] . 

Having derived the general expression for transmission 
of electromagnetic fields through a subwavelength chan- 
nel containing, to this point, conventional material pa- 
rameters, we now examine the effects of taking the limit 
of vanishing e. It is relatively straightforward to show 
that Eq. I|10p for the electric field at the slit exit now 
reduces to the succinct form, 
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where H z = 2qIoE^ ot , and is spatially constant in the 
slit. It can be deduced from Faraday's Law, H = 
— i/(/ifco)V x E, that E x must therefore be a linear 
function of the coordinate y within the slit: E x = 
(y/l)SE x + E x °p 7 where SE X = E x °v - E x ot . The elec- 
tric field at the top of the slit, £^ op , is simply, £*°p 

1 - 2iqI k ln)E^ ot . For q < 1, we can write Eq. {Igj) 
strictly in terms of elementary functions, 



2ir[ir + k lnq(3 -2^-2 big)] 
4 ln<?(27 - 3 + In q) + tt 2 + (2 7 - 3) 2 ' 
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demonstrating that for a given ratio of slit width to wave- 
length, t is simply a linear function of the channel length 
I and permeability (i. 

To illustrate how the transmission depends upon the 
geometry of the slit, we present in Fig. [2] a three- 
dimensional representation of r [Eq. (|13|) ] as a function 




FIG. 2: Transmission, r, through a zero-e slit as a function 
of the dimensionless geometrical parameters, l/X and p/A. 
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FIG. 3: Transmission characteristics for a zero-e metamaterial 
slit, demonstrating the / and p that give the optimal transmis- 
sion at a given wavelength. Brighter regions indicate higher 
transmission. Relatively high transmission can persevere in 
deep slits (I S> A) and small openings (p<A). 



of the slit dimensions, I and p, scaled by the wavelength. 
In this figure, the source field is normally incident with 
wavelength, A = 5/j,m, which is very close to the res- 
onant plasma wavelength, so that the frequency disper- 
sive e nearly vanishes. It is evident from the plot that the 
transmission does not possess Fabry-Pcrot-like resonant 
oscillations as a function of I, as would be expected from 
a slit with conventional material. The figure shows that 
energy flow is generally restricted unless the geometri- 
cal parameters coincide with the brightly peaked regions 
that can decay quite rapidly. 

The wave impedance, defined as Z = E x /H z , is cal- 
culated from the field expressions above to yield the fol- 
lowing impedance relation between the top and bottom of 
the openings, Z top = Z hot -ik l^, with Z hot = l/(2ql ). 
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FIG. 4: Transmission through deep metamaterial slits as a 
function of p/X. A considerable resonance peak emerges for a 
small subwavelength width of zero-e material. Also shown is 
a matched zero index medium where p, and e are both nearly 
zero. The inset shows the intersections of the lines with the 
y(q) curve, corresponding to electrical widths that yield peak 
transmission. Differing ratios of l/X (see legend) yield lines 
with different slopes. These results are consistent with the 
main plot, where the channel length corresponds to I = 50 A. 

A remarkable property of a narrow zero-e slit is that de- 
spite the considerable wave impedance for p <C A, en- 
ergy can still be transmitted if the slit length I is in- 
creased. This behavior is further illustrated in Fig. [31 
where the transmission is shown to be robust for an ex- 
tended range of channel lengths and widths. The maxi- 
mum transmission, r max , seen contained within the con- 
tinuous bright curve, is determined by examining the ex- 
trema of Eq. JI3]). By taking the appropriate deriva- 
tive of the denominator, we find that r max is determined 
by the transcendental equation, fiqkol — 7(g), where 
■y(q) = — Im[/ ]/|^o| 2 is a gradually increasing monotonic 
function of q. Thus, as a function of q, and for each 
the intercept of "f(q) with lines of slope /ih^l yields the 
permissible electrical slit widths that achieve peak light 
transmission. This technique is shown graphically in the 
inset of Fig. [JJ where the opposite relationship between 
slit width and length are clearly seen. 



It is also of interest to determine if the transmission 
can be enhanced for light incident upon a a matched 
zero index material, whereby both e and /1 are vanish- 
ingly small. In the main plot of Fig. [¥J the transmission 
is therefore plotted as a function of slit width for a zero-e 
slit (with \x = 1 as usual), and matched zero index slit. 
Surprisingly, for zero-e media, the shown subwavelength 
widths support transmission resonances that are absent 
in matched zero index slits. Although the I dependence 
washes out for matched zero index channels, r still de- 
pends strongly on the width dimension, and Eq. (|13p ap- 
proximately reduces to the expression for a slit filled with 
a conventional dielectric and length satisfying the Fabry- 
Perot geometric resonance condition. We can thus con- 
clude that for zero-e media and narrow enough openings, 
energy flow can be much greater than for slits loaded with 
metamaterials having no wave impedance mismatch be- 
tween the ends of the slit. We also considered the effects 
of finite T, and concluded that the presence of absorption 
loss simply reduces the overall magnitude of r by a fac- 
tor that becomes greater with increasing slit depth, but 
having no effect on the geometrical parameters leading to 
the particular resonance location. This bodes well for fu- 
ture development of miniaturized devices used for energy 
transport and narrow band frequency selective surfaces. 

In conclusion, we have shown that significant transmis- 
sion arises for a range of subwavelength widths in suffi- 
ciently deep slits containing zero-e inclusions. We also 
demonstrated that these resonant channels can transmit 
greater energy compared to the corresponding matched 
zero index slits. Although there is currently limited 
metamaterial fabrication in the IR and optical frequen- 
cies, some recent progress has been made with negative- 
index metamaterials involving metal layers separated by 
a dielectric [l7[ , alternating layers of InGaAs and AlInAs 
semiconductor s llal . and photonic crystals in the near- 
infrared range [l9(. 
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